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ABSTRACT 


"  A  general  variational  approach  to  study  systems  composed  of 
complex  charged  molecules  is  discussed.  In  this  approach  the  variational 
trial  functions  for  the  free  energy  functionals  are  constructed  from  the 
asymptotic  limiting  (AL)  forms  of  the  direct  correlation  functions. A 
number  of  examples  are  discussed,  and  in  each  case  the  variational  form 
of  the  direct  correlation  is  given  explicitely.  The  relation  to 
Onsager's  procedure  of  immersing  the  system  in  an  infinite  conducting 
fluid  for  obtaining  an  energy  bound  is  discussed  in  detail. 


:v«v '-f 


Introduction 

Real  fluids  are  composed  of  molecules  that  are  objects  of 
complex  geometries  and  charge  distributions.  In  a  previous 
note^^  we  have  shown  that,  by  studying  the  asymptotic  high 
density  limit  ( AHDL)  and  the  asymptotic  strong  coupling  limit 
(ASCL)  one  is  able  to  reduce  the  problem  of  computing  the 
thermodynamics  and  correlation  functions  of  the  system  to  a 
geometrical  calculation  involving  overlap  integrals  between  the 
objects. 

'In  previous  work  *  a  simple  geometrical,  physically 
intuitive  meaning  of  the  direct  correlation  functions  (dcf)  for 
point  charges  in  a  background^ (as  interactions  between  smeared 
charges)  and  hard  spheres  (as  overlap  volumes)  within  the  mean 
spherical  approximation  (MSA)  was  given,  thu3  also  revealing  its 
analytic  structure.  As  a  result,  the  above  program  can  be 
carried  out  completely  for  relatively  simple  systems  (as  e.g.  the 
general  ionic  mixture  of  the  multicomponent  plasmas^’^)  using  the 
MSA  free  energy  functional  which  interpolates  between  the  exact 
weak-(Debya-Huckel)  and  strong-( "Onsager-ty pe" )  coupling  bounds 
for  the  potential  energy.  Though  featuring  fewer  "idealistic#" 
features,  in  view  of  the  higher  complexity  of  the  problem,  this 
approach  was  succesfuly  used  to  analyze  the  "isotropic"-"nematic" 
transition  of  line-charges^ ’  and  the  coupling  of  the  growth  of 


micelles  to  their  degree  of  alignment0.  In  the  present 


3 


communication  we  extend  these  methods  to  a  much  larger  class  of 
objects . 

The  proposed  approach  is  to  write  down  an  approximate  free 

energy  functional  which  has  to  be  variational  with  respect  to  the 

pair  functions.  These  would  be  either  the  indirect7 
8  9 

or  direct  *  (i^jCr^))  correlation  functions.  In  order  to  get  a 

convenient  formalism,  we  have  to  use  simple  functions  with 
physically  motivated  coefficients.  Inde*^  the  direct  correlation 
function  in  the  asymptotic  limits  ( AL  ■  either  AHDL,  ASCL) 
provides  such  a  simple  intfructive  basis. 

The  approximate  solutions  also  provide  exact  bounds  for  the 
free  energy  of  system.  In  the  present  work  we  present  a  few 
results  for  simple  system. 

In  section  2  we  discuss  the  Mean  Spherical  Approximation 
(MSA)  for  hard  core  -Green  function  systems.  The  ASCL  of  the  MSA 
and  Hyphernetted  Chain  (HNC)  ar e . d i sc u ssed  in  section  3.  In 
section  4  we  give  the  general  "Onsager"  solution  for  the  charge 
smearing  problem,  which  provides  the  basis  -  set  functions  for 
the  MSA-HNC-ASCL  problem.  The  variational  solution  by  expansion 
in  this  set  is  discussed  in  section  5.  Section  6  is  devoted  to 
the  discussion  of  the  PY  theory  for  hard  objects,  in  parti;ular, 
the  scaling  of  the  direct  correlation  function.  An  interesting 
application  to  bonding  and  aggregation  within  the  MSA  is  given  in 
section  7.  Section  8  contains  a  general  discussion  and 
cone lusions . 

c 

r 

i 

r. 

r. 
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2.  Mean  Spherical  Approximation  for  hard-core.  Green  s  function 


systems 


The  general  charged-hard-ob jects  system  to  be  considered 
(and  termed  hard-core-Green ' s-f unction  ( HCGF)  system)  consists  of 
hard  objects  with  imbedded  charge  distributions 


►  ->  . 

p  <•:) 

N. 


q.  f  <■£' 

^i<s  1  - 1  & 


(2.1) 


Ou  * 

where  p.  .  -  isV^ormalized  charge  distribution  of 

J  v  A  G 

multipolarity  ^  ,  associated  with  the  Green's  function  (GF)  G, 

and  o^.^s  corresponding  coupling  constant  (•"charge").  In 

this  discrete  representation  for  polydispersity ,  an  object  i  of 
relative  concentration  X^is  considered  distinct  fron  j  if  one  of 
the  characterizations  in  the  tuple  i  ■  (shape,  size,  orientation, 
charge  distribution)  is  different  from  that  in  j.  In  addition  to 
the  hard-core  repulsion,  the  pair  interaction  between  two  objects 


i,  J  at  a  distance  £  is 


4>..cO  *  II 

7vi  ua. 


e<s,  e  g, 


(2.2) 


where 


( 

For  example,  a  point  monopole  and  a  point  dipole  at  the  center  of 
an  object  are  represented  by  p.  =  q,.  c> 


and  C  respectively  .  The  Coulomb  and 

screened-Coulorab  potentials  (  *  ‘/i*i  ^  ^  «•*  eLr  'V*1  VlM 

respectively  in  three  dimensions)  are  the  GF  for  the 

and  ^  Liouville  operators.  We  use  the  following 

d 


notations: 


-  >  total  number  density  ;  /3-»  '/v*  ' 


inverse  temperature;  overhead  tildes  denote  Fourier  transforms; 

"det"  denotes  the  determinant; 

3)  ■  dimensionality;  -  unit  matrix  dl.  ;  C  -  matrix  of 
dcf 's  with  elements  *->  «  51  *•  «. 

‘i  '  'I*  “  k  Ko  (S 

total  monopolar,  G  -  type,  Charge  density;  m  0,1  according 

*T  ** 

the  whether  (k’-o)  is  finite  (e.g.  Yukawa)  or  infinite  (e.g. 

Coulomb) . 

The  MSA  equations  are  obtained  from  the  Ornstein-Zer nike 
(0Z)  relations  between  the  direct  correlation  functions  (dcf's)  c  C 
and  pair  correlation  functions  (pcf's)  e» 

wi  J 

W.  (t)  *  ,  w  ^  x  J-')  S/-'1 

'•i  »  (2.4) 


and  the  closure  relations: 


a  (l\  z  O  r  C  <S[- 


(2.5) 


for  the  inner  hard  core  excluded  region,  and 


V-'  I7  5>; 


(2.6) 


J  J  J  / 


core 


Let  denote  the  value  of  f~  for  B-0,  i.e.  for  the  hard 

O  ' 

system. Recall  that  represents  the  compressibility  factor, 


Z-  pV/NkT 


as  obtained  from  the  PY  (i.e.  MSA)  equation  for  the  hard  core  system, 


via  the  relation  [9,10] 


F  -d/2)  ( 2  -1) 


(2.10) 


where  Z  denotes  Z  as  obtained  from  the  compressibility  equation  of 


state: 


K>p/4nyj*'£'c-(~ 


(2.11) 


If  f  represents  the  excess  free  energy  per  particle  in  units  oi 
kT  and  f  =■  lim  f  ,  then  the  MSA  approximation  states  that  [10] 


f  -f  -  F  - 


(2.12) 


I 


for  the  equation  of  state  obtained  from  the  expression  for  the  energy 
the  system  : 


u 


/ 


-U/NkT-l/2  n 


g  (r) 


(2.13) 


Central  to  our  treatment  below  is  the  "Ewald"  identity  for  an 

function  (Ewald  function)  9..(r)  for  which  the  Fourier  transform  ©  (fe) 

h  - 

exists  [4,10]  : 


1/2 


n  }  x.x.  dr  g  (r)  9.  (r)  ■ 

Ch  k  j  J  ~  ' 

S  1/2 


n  ;  x.x 


dr  9  (r) 


-  1/2  (2\vj^V  jdk  9*  (k)  + 


1/2  (in^Rx.  xj1 


dk  9.  (k)  S.  (k) 

~  “  S  * 


(2.14) 


where  the  structure  factors  S..(k)  defined  by 


(2.18) 


(2.19) 


u»  B[C]  + 1/2 


fom  which  (2.11)  is  obtained  for  systems  of  hard  core  objects  upon 

setting  u  =0. 

*s»» 

Finally  note  that  when  the  total  monopolar  charge  of  the  system 

/•w 

is  not  zero  we  need  ,  in  the  case  when  <f>  (k*0)»0  (e.g.the  Coulomb  case 

<* 

to  introduce  a  compensating  background  charge  density  which  is  uniform 
i m  all  space,  i.e.  it  penetrates  the  hard  objects.The  background 
density  is  equal  to  -  p  .  The  potential  energy  of  the  system 
including  the  background  is  given  by 


U/NkT»l/2 


n  x-  k  Jdr  h.  (r)  B  & 


(r) 


(2.20) 


The  general  expression  for 


is  thus  given  by 


_p^p/^i-  n  jd~  [C4i(E> 


while  the  functional  B  is  written  as 


(2.21) 


The  sum  over  orientations,  implicit  in^x  x.  ensures  that  for  th 
HCGF  systems  * 

Zx  x  (dr  fa  (r)  »0  (2.23) 

C  J  J  ~  I'i  * 

whetnever  total  monopolar  charge  neutrality  is  preserved. In  view  of  thi 
(see  3.11)  below),  (2 . 20) - (2 . 22)  contain  (2.8), (2.13)  and  (2.17)  as 
special  cases. 


a)  Diagonal i2at ion  of  the  variational  free  energy 


A  strong  coupling  (SC,  superscript  oO  )  limit  for  the  MSA 

is  reached  when  either  the  hard-core  compressibility  tends  to 

zero,  or  when  any  of  the  coupling  constants  (the  charges)  tend  to 

infinity  so  that  the  free  energy  is  dominant  by  the  energy  term, 
P-*  -o 

j  ■  B.  In  either  case  the  dcf's  diverge  ,  and  in  order  to 
satisfy  the  non-negativity  of  the  argument  of  the  logarithm  in  L, 
the  diverging  dcf's  satisfy  ' 


r  c'wf  =  ^  c.- 

^  v.  j  i 


c“  it) 


(3.1) 


(3.2) 


In  view  of  the  closure  (2.6)  which  servers  as  a  boundary 
condition  for  the  GF-potentisl  (2.2;  2.3)  these  relations  can  be 
satisfied  only  by  the  following  convolution  type  forms 


constant  •  P  P  ‘(.V)for  hard  core  SC  (3.3) 

S~,tl  ■ 


•~0i\  -  f  °0  ^ 

V1  '  Ji~". 


for  charge  SC  (3.4) 


where  SC  denotes  strong  coupling. 


The  normalized  "smearing"  d  i  s  t  r  i  l>u  t  i  o  ■■■  •»  .  O  t  A  '  o')  •=.  \ 
are  confined  to  the  volume  of  the  hard  object  i  or  its  surface 
and  must  satisfy  the  MSA  "boundary  condition"  (4).  For  (3.4)  we 
specifically  consider  the  new  ("smeared"  as  opposed  to  smearing) 


v 


i  v 


coupling  constant  .  such  -  n  a  w  w 

l  (j| 

interaction  between  objects:  and  J 


<<5,  C.V 


<-0 


j 


(3.5) 


will  satisfy  the  MSA  boundary  condition 


(  \  *  -  x1  -  rij 

e<s .  e‘<s'  ~ 


(2.6),  n  fil  y 


Z  (s.  Z  G 


"  ti 


Z.C,,*  Cr 


r  7  <r.  . 

(3.6) 


Using  (3.3)  and  (3.4)  we  now  consider  the  strong  coupling 
MSA  problem,  which  may  be  posed  separately  for  each  diverging 
component-  i.e.  either  hard  core  SC  or  the  SC  limit  of  any  of  the 
independent  coupling  constants 

90  y$ 

Inserting  (3.3)  into  B  r3  [  (2  ]sO  we  find  a  du-^ct 
generalization  of  the  hurd-spheres  result^  ,  namely  that  the  MSA 
hard-core  EOS  diverges  when  the  total  packing  fraction, 

^  ■  (volume  of  objects/volume  of  system)  ■  1  with  the  dcf's 

satisfying 

Cf.  10/4*° At'-')  ^.V  tth 

^  ^  (3.7) 


overlap  volume  of  objects  i  and  j  with  separation  r_  divided  by 
the  overlap  veluevfiat  zero  separation. 


In  view  of  the  additivity  of  the  total  interaction  potential 


in  terms  of  the  independent  "charges",  we  consider  at  once  the 
most  general  case  of  strong  coupling  of  all  charges,  although  as 
stated  above,  each  may  be  considered  separately. 

Inserting  (3.4)  into  the  SC-MSA  functional  B  [ C *°  ]  we 
obtain  the  Onsager-ty pe^ ’ ^  expression  for  the  lower  bound  for 
the  true  potential  energy  of  the  system,  given  by  (for  reasons 
given  below  we  prefer  the  "smeared"  new  charges  g  o xJ^-r 

This  is  choice  of  convenience  and  using  either  gives  identical 
results). 


y  ft  L  c"  1  ~ 

P 


,C  Q.  Q 


(3'8) 


2. 


CC'-o') 
S.&,  £ G 


•.  *e 


(3.9) 


The  first  (monopolar)  term  in  (V^),  an  Onsager  background  term, 
vanishes  for  globally  neutral  systems  (  o  -  O  ).  The 
second  terra  is  minus  the  self  energy  of  the  "smeared"  charge 


distribution.  Note  the  decoupling  o£  all  the  components  in 
(3.8),  which  is  typica i  ~  to  the  Onsa&cr  best  bound  scheme  and 
is  analogous  to  the  diagonalization  of  the  Hamiltonian  (sec  Sec 
v) .  The  new  "diagonal"  MSA  problem 


<S  <? 

Ji  e  & 


(3.10) 


becomes  an  electrostatics  problem  o  f  finding  the  optimal  On  s  a  g e  r 

.smeared"  distributions . ^he  derivation  of  (  3  .  &  1  ,  (  3  .  9 )  is  made 
transparent  in  (b)  below. 


Even  though  a  comprehensive  analysis  of  the  HNC  theory  for 
hard  objects  is  the  asymptotic  strong  coupling  limit  (ASCL)  has 
not  been  performed  yet,  the  charge  strong  coupling  limit  of  the 
HNC  can  be  shown  to  be  identical  to  that  of  the  MSA.  Thus,  the 
results  (3.  8  )  and  (3.10)  are  equally  valid  for  both  HNC  and  MSA 
and  the  corresponding  "best  bound"  problem,  which  as  we  show  in 
the  next  section  has  been  solved  already  by  Onsager,  many  years 
ago,  provides  the  ASCL  for  charge  strong  coupling  for  both  HNC 
and  MSA. 


(b)  The  Evald  identity  and  the  Onsager  Process 


Considering  the  Evald  identity  (2.14),  notice  that  electrostatic 

interactions  between  charged  particles,  namely  d  (r)  or  W  (r)  of  (3.6), 

'  iy  ^  » *-i  - 

are  legitimate  Evald  functions  9.  (r) .  The  left  hand  side  (l.h.s.) 

‘•j  *■' 

of  (2.14) is  correpondingly  the  total  pair  interaction  potential  energy 
per  particle,  of  the  objects  carrying  the  original  charge  distributions 

O 

p^(*)  and  the  smeared  charge  dstributions  p.  (x) .  The  first  term  on  the 
right  hand  side  (r.h.s.)  of  (2.14)  vanishes  .for  an  electrically  neutral 
system  (  Tx.  q.  =*0  )  ,  since 


y  x  x.  (dr  (r)  ■  lim  Vx.x.q.q  p°  (k)  p°  (k)  ^  (k) 


> 

*  lim  (Tx.q.p*(k)  )  i  (k)  (3.11) 

v 

and  similarly  for  ^x.  x  jdj  (r)  .  (k)  is  the  FT  of  the  GF  potential 

(e.g.  Coulomb)  .The ’‘second  term  of  the  r.h.s.  of  (2. 14)  is  easily 
recognized  to  be  minus  the  self  energy  of  the  charge  distributions,  so 
that  the  last  term  of  the  r.h.s.  of  (2.14)  represents  the  total 
electrostatic  energy  (per  particle)  of  the  system  of  charged  objects: 


(3.12) 


whereuj  >o  is  the  surface  of  a  D-dimensional  unit  sphere,  E  and  E 


are 


V  ->*»  ,  N/V-»n)  . 

The  Onsager  process  [11]  is  equivalent  to  adding  and  substractir.g 
the  total  electrostatic  energy  of  t.*e  smeared  charge  distributions,  and 
rearranging  terms. Using  (2.14),  it  is  equivalent  to  the  following  Ewald 
identity  (see  Fig.  1  ): 

U/N-  (  U/N-  U  /N)  +  U  /h‘ 

(3.13a) 

i.  e. 


For  hard  core  Green's  function  systems  (HCGF)  (2.5)  is  exact,  so 
that  due  to  (3.6)  the  first  term  on  the  r.h.s.  of  (3.13b)  vanishes. 3y 
(3.12)the  second  term  of  the  r.h.s.  of  (3.13b)  is  non-negative,  so  that 
the  Onsager  type  bound  is  obtained: 


U/N  > 


(3.14) 


where  u.  is  the  self  energy  of  the  smeared  charge  distribution  q  p(.x) 

i)i 


arbitrary  N  (nur.be r  cf  cb]ecrs)  provided  than  one  elecoror.euora . . o  .■ 

preserved. 

Returning  to  the  SC-MSA  problem,  and  denoting  by  "overbars"  the 
optimized  quantities  obtained  from  the  solutions  of  (3.10),  Ve  finally 
get: 

i)  The  SC-MSA  result  is  an  exact  lower  bound  to  the  potential 
energy  of  the  system. 


U/N>  (  U/N  ) 


MSA  free  energy 


■  the  interaction  between  the  optimally  smeared  charges  in  objects 
i  and  j  of  separation  r  . 

(3.16) 

Note  that  as  stated  in  "words",  namely  in  terms  of  the  basic 
characteristics  of  the  interactions,  our  simple  expressions  (3.15)  and 
(3.16),  uncover  th«  Mf5A  meaning  for  arbitrary  dimensionality. 

When  there  is  a  background  (3.14)  takes  the  form 


-  l/p  B[C  “°  ] 


(3.15) 


Onsager  Bound, despite  the  approximate  nature  of  the 


Note  however  that  the  second  term  of  the  r.h.s.,  which  represents 
the  difference  of  the  interaction  energies  of  the  smeared  charges  and 
the  original  charges  with  the  uniform  background,  is  also 
diagonal .  Indeed,  we  may  write  ((featuring,  in  full  i<*->i  ,  £,G) 


1/2  n  2.  /dr  [ip.  (r)-£.  (r)  ]- . 


•1/2 


n]x;  x<5.  <3  lim  CP  „(*)p  (k)-  p*  (k)  P  (k)  (k) 

L,  J  l'<s  J06  -  >4  -  * 


£ . 

to 


ce.i. 


■lira  1  /2  n  ()x.  q  [p  (k)+S°  (k)])  (Tx.q  [p  (k)  -  pJ  (k)]d(k)  £ 

^  L,  *•  '•oG  J,;  ~  )  ”  <-,J  Jsi  VS  rj  o  6  "  'Q  iJ 


>■ « 


1/2  fr4Ixil/qi^{d2 


(3.18) 


and  we  recall  that  p(k-0)-p*  (k-O)-l. 


Comparison  with  (3.9)  shows  that  the  MSA  or  HNC  liquid  theories, 
dictate  that  the  smeared  system,  represented  by  U  /N  in  (3.13a)  and 


co 


discarded  in  obtaining  U  /N  ,  should  be  charge  neutral  also  when 

— * 

k"°)  .That  is  the  reason  for  ^  in  (3.9)  when  ^  (k-0)-  o  and  we 
formally  do  not  need  to  introduce  the  uniform  background,  since  the 


Or.sager  process,  dictates  that  the  "s.-.eared  syster."  snculd  cc.tta.r. 

Uniform  Backgroundknd  thus  be  totally  charge  neutral.  Since  U  i 
discarded  in  the  MSA-Onsager  estimate  cf  the  potential  energy  of  th 
given  system,  it  can  be  expected  to  be  relativley  small  only  for  a 
totally  charge  neutral  system. 


4 . "Onsager  Solution"  of  the  MSA  or  HNC  for  Charge  Strong  Coupling 

In  his  classic  paper  of  1939,  Onsager  [11]  considered  the  problem 
of  obtaining  a  lower  bound  to  the  potential  energy  of  a  system  of 
charged  hard  objects.  His  method  of  solution  is  consistent  with  the  MSA 
boundary  condition  (3.6)  and,  in  fact,  provides  the  solution  of  the  ASCL 
for  both  the  HNC  and  MSA  (3.10) .We  consider  in  this  section  the  Coulomb 
potential, to  avoid  the  complications  and  specific*  details  of  other  gf 
potentials,  to  which  the  treatment  also  applies. This  will  allow  us  to 
use  elementary  electrostatics  in  3  dimensions. We  discuss 
electrically  neutral  systems,  ys  -0,  in  the  first  place,  then  systems 
that  are  not  electrically  neutral,  and  therefore  require  a  neutralizing 
background  for  thermodynamic  stability  and  at  the  end  of  this  section  we 
give  examples  for  the  special  case  of  centrally  charged  hard  spheres. 

a)  Total  Charge  Neutrality  , 

There  are,  in  general,  an  infinfite  number  of  ways  to  replace  the 

charge  distribution  of  an  object , p  (x)  by  a  smeared  distribution 

with  the  same  potential  outside  of  the  object.  The  most  obvious  ones  are 

those  associated  with  the  spherically  symmetric  distribution  p  (|xi) 

1 1  «1 » <  (, 

which  by  Gauss's  theorem  satisfy  the  required  boundary  condition 
(3. 6). There  is  ,  however  only  one  surface  distribution  jo.(^(x) 
which  satisfies  (3.6).  This  surface  distribution  is  the  one 
corresponding  to  minus  the  induced  charge  on  the  surface  of  the  objects 
when  the  system  is  immersed  in  a  uniform  conducting  fluid,  as  originally 
proposed  by  Onsager.  As  is  discussed  below  in  a  special  case,  this 


From  this  solution  we  compute  the  component:  of  me  elec:::: 
normal  to  the  surface  of  the  object,  and  just  inside  that  surface 


.  (x) .n| 


(4.7) 


^  :  S 


The  induced  charge  on  the  surface  of  the  object  is  given  by 


S'  (s)—  E.\(s)/4 

V.  ** 


(T 


(4.8) 


The  Onsager  solution  for  the  optimal  smeared  charge  distribution 
J.  (x)  is 


Onsager 's  original  bound  to  the  potential  energy  of  the  system 
was  given  [11]  in  terras  of  the  total  potential  energy  of  the  "Onsager 
objects" , namely  the  neutral  objects  consisting  of  the  original  charges 

O 

p.  (x)  plus  the  surface  induced  charges  induced  by  them,  O’,  (s)  .  It  is 

y  h  *•  ^  /v 

easy  to  show  that  this  energy  is  equal  to  minus  the  self  energy  of  the 
surface  charge  distributions^ 


I  A  (s)  O'(S)  +1/2  (d 

s  (s)  o'(s) 

(4.10) 

J 

'  So  -  ~  J 

~  v 

from  ( 4 . 4 ) , however  we  obtain 


■  ids  p  (s)  o(s)  =  -(self  er.eray  of  the  surface  charoe 

J  -  -  - 

distribution  d(s))  (4.11) 


In  order  to  get  some  feeling  for  the  optimization  of  p. (x) 

'  W  >* 

obtained  by  the  surface  distribution , consider  the  special  case 
i  (s)=»c  “constant  (e.g)  a  point  monopole  at  the  center  of  a  SDhere)  . 

The  interior  potential  is  given  by  (x)  =  ^(xj-c^,  i.e. 

<^(x)  Cc,  the^lectri c  field  inside  the  object  due  to  the  surface  charge 

is  zero,  EC=«O.OutsideJthe  objective  still  have  from  (4.6)  that 

o  •> 

E.  =E.  .The  general  conditionfor  the  validity  of  (3.6) is  that  the  field 

*»V(0 

outside  the  object  due  to  the  smeared  charge  is  the • same  as  that  due  to 
the  original  charge, 


?  \  ^ 

.  ]  -E  . 

v  ) 

e*r  «4 


(4.12) 


The  self  energy  of  the  smeared  charge  distribution  i  given  by 


U.  -1/S  tfjr-oj-i/airJdVKE.^^^i1 

ati  )|M( 

“1/8 K  (dV|(E<)  f  +1/87  fdV|  (E  )  |"* 

J  ‘  **«c**«A  J  v*.c 


o'o^c.t 


outuM* 
o'o  j 


(4.13) 


The  absolute  minimum  for  this  quantity  is  reached  whenfE4^ 

$  *  t 

>  which  in  the  case  considered  above  is  obtained  by  the  Onsager 

solution  with  E^*0  for  £  (s)»  constant. 

ri  •  - 


The  Onsager  procedure  for  the  MSA-HNC  best  bound  problem  of 


immersing  the  system  in  a  uniform  conducting  fluid, works  well  when  t 
is  a  neutralizing  background, The  Or-_ager  object  now  consists  of  the 
original  charge  distribution  jO°(x),  the  uniform  background  density 
"trapped"  inside  the  object  ,jo  =-ja  /  and  the  induced  charge  on  the 
surface  due  both  to  ^  and  .Eg.  (4.1)  now  takes  the  form 


t.  ,<• 


\?  f  (x) 

1  v  ** 


with  the  condition  (4.2).  Denoting  by  ^  (x)  the  solution  of 


(4. 


V  ^  (x)  *“4\\o  inside  object 

h'0  b 


V  f.  (x)  =0  outside  object  (4. 


with  the 
(4.3)  to 


condition  that  4  (x-*«©)-0,  we  use  <&  (x)  as  defined  prior  to 
solve  (4.3)  with  the  boundary  condition 


to  get 


(4.16) 


(4.17) 


Thus  4  (xV  ,4  (X\  and  4  (X)  are,  respectively,  the  electrostati 
/ i,o  -v. 

potentials  due  to  d0)  ,p  and  O*  (s)  in  the  Onsager  object.  Let  us 
denote  by  o  the  original  fixed  charge  of  the  object, by  b  the  backgrour. 
charge  and  by  s  the  surface  charge,  and  ob,os,bs,bb, ss, . . the 
electrostatic  energies  between  those  charge  distributions  .bb  and  ss  a 
the  self  energies  of  the  distributions  jo^  and  0"  (s).It  is  easy  to  see 
that  (3.18)  may  be  written  ,  in  the  case  of  Onsager  smearing  as 


2 


x . (bo+2bb+bs) 


(4.19) 


while  the  first  terra  in  (3.17),  the  self  energy  terra  of  the  components 
in  ®  .  , is 

1  «•  V 


.  . (r-0) 
*-  C 


•  (ss+  bb+bs) . 
v.  \. 


V 


(4.20) 


thus  ,  (3.17)  takes  the  form  (see  Appendix  B) 
cO 

U/N  > (U/N)  - 


x-  (ss-bb-bo). 


(4.21) 


Consider,  in  turn  the  total  energy  of  the  Qnsager  object 
(relative,  as  usual  to  the  self  energy  of  p° ) : 


u  * (bs+bb+ss+bo+os). 


(4.22) 


It  easy  to  see  from  (4.16)  that 


bs+os«-2ss 


(4.23) 


so  that  (4.21)  becomes  identical  to 


£.(U/N)  xu. 

L  1 


*«•  *■ 


(4.24) 


Thus,  when  there  is  a  background,  the  Onsager  solution  for  the 
optimal  smeared  charges  is 


f;  <?>  *-  <'  <;>  -fb 


(4.25) 


(c)  Special  Case:  hard  spherical  centrally  charged  objects. 
The  most  general  system  for  which  a  formal  solution  to  the  MSA 
has  been  found  consist  of  hard  spheres  of  arbitrary  sizes  with 
central  point  multipolar  charges,  for  which  the  solution 

to  the  MSA-HNC-ASCL  problem  (3.10)  can  be  summarized  by 
the  following  simple  statement  (see  Appendix  C): 


All  multipolar  charges,  [  )  1,  should  be  uniformly  smeared  on 

the  surface  of  the  sphere.  The  monopolar  charge,  a.  .should 

U  o  (3 

be  smeared  uniformly  on  the  surface  and/or  uniformly  in  the 
volume  of  the  sphere  according  to  whether  ^'->3  or  j3.  &  f  0 

respectively . 

All  the  analytic  solutions  of  the  MSA  for  HCGF  systems,  that 

12  13 

are  known  to  us,  including  hard  spheres  *  ,  hard  spheres  with 

1 4 

centrally  imbedded  point  monopolar  charge  (with  and  without 

background ^ f  hard  spheres  with  centrally  imbedded  dipoles^  and 

17  IS 

multipoles  ,  mixtures  of  charged  hard  spheres  with  point  ions  , 

19 

hard  spheres  with  monopolar  Yukawa  charges  ,  obey  our  relations 
(}•  ~1)  1  (i.  f  0  and  ^.\du).  with  the  analytic  forms  given  by  (  J.l  l\(b.  12) 
maintaining  it  at  all  densities  and  temperatures  (couplings). 


5.  Variational  Solutions 


of  the  MSA  by  Expansion  in  the  Onsager 


"Basis  Functions". 

(a)  Analogy  with  the  Hamiltonian  variational  problem. 

A  typical  liquid  state  calculation  involves  a  (usually 
approximate)  free  energy  functional  which  has  to  be  variational 
with  respect  to  the  pair  functions.  In  analogy  to  the 
Hamiltonian  (  free  energy)  ground  state  energy  problem,  the 

following  "ideal"  situation  is  desirable: 

(i)  To  have  a  physically  motivated,  good  choice,  of  a  set  of 
basis  "wave"  functions  (  9  pair  functions)  which  obey  all  the 

symmetries  of  the  "hamiltonian"  (  <f — =?  geometry  of  the 

object  and  the  charge  distributions):  so  that  (ii)  the  pair 
functions  could  be  expanded  in  terms  of  the  basis  set,  with  only 
the  coefficients  (which  can  be  assigned  intuitive  physical 
meaning)  to  be  determine^by  the  variational  free  energy.  A 
mlnimun  set  of  coefficients  must  still  involve  the  geometry  of 
the  objects  and  the  values  of  the  charges  (iii).  Eventually 
increasing  the  size  of  the  basis  set  (the  number  of 

coefficients),  -tore  accurate  solutions  can  be  obtained. The  exact 
solution  is  reached  when  the  complete  set  of  basis  functions  is 
involved. In  general,  however, a  small  number  of  coefficients  will 
be  required. 

( 1  v ) A,  la  th.  Haniltonlaa  problea  it  la  daalrabla  to  obtai, 

fr°»  Mfr.,l..tj  expansion  and  ex*=t  bound  to  th, 
the 


system. 


energy  0f 


The  general  idea  of  solving  che  variational  MSA  problem  by 
using  a  trial  dcf  with  coefficients  to  be  determined  by  the 
variation  equations  in  certainly  not  new.  Ne’er  before  we  had, 
however,  at  our  disposal  the  complete  set  of  basis  functions  as 
we  now  have  due  to  the  mapping  of  the  ASCL-MSA  on  the  Onsager 
scheme,  which  provides  the  exact  analytic  form  of  the  MSA  dcf's. 
With  trial  dcf's  having  the  exact  analytic  form  of  the  full  MSA 
solution  it  is  only  a  method  of  employing  the  full  set  of 
required  coefficients  (i.e.  expanding  the  solution  in  the 
complete  set)  in  order  to  have,  from  the  variational  equations 
for  the  coefficients,  the  exact  solution  of  the  MSA.  The  Onsager 
basis  set  of  functions,  constructed  by  (3.12)  using  the  Onsager 
smeared  charges,  (i)  h<ts  a  physicaly  intuitive  meaning,  (ii) 
corresponds  to  an  exact  lower  bound  to  the  potential  energy,  and 
(ill)  is  constructed  by  using  elementary  electrostatics 

employing  the  basic  geometrical-physical  constraints  that  must  be 
taken  into  account  -  the  liquid  state  theory  part  of  the  problem 
does  not  make  the  problem  more  complex  than  it  already  is  at  the 
basic  electrostatic  level.  Even  if  we  "stick"  to  the  full  basis 
set,  and  attempt  an  exact  solution,  this  procedure  has  the  clear 
advantage  of  providing  a  direct  physical  description  of  the 
solution . 

We  may,  however,  with  increasing  experience,  use  only  those 
elements  of  the  basis  set  which  are  more  important.  The 
conceptual  analogy  with  the  widely  used  Hamiltonian  variational 
solutions  is  complete. 


(b)  Practical  details:  Examples  and  Soraments. 

The  complete  basis  set  of  functions  for  the  MSA  consist  of  the 
overlap-volume  functions  (3.7)  and  the  Onsager-smeared 
interactions  (3.1.6).  Recall  that  the  MSA  dcf's  outsiae  the  cores 
are  already  specified  by  the  closures  (2.6),  so  that  without 
further  mention  it  is  understood  that  only  the  region  inside  the 
core,  enters  into  the  discussion. 

Exaraple(bl):  restricted  primitive  mo<iet(RPM). 

The  RPM  consists  of  a  binary  mixture  of  equal  sila  and  opposite1.-; 
charged  spheres.  The  overlaps  volume  function  is  of  the  form  Lil 


Chc(x)  "  Ao  +  Alx  +  A3*3 *  x  *  1  <  1  ^5.1) 

The  Onsager-smeared  interactions  between  two  uniformly  surface 
charged  spheres  is  of  the  forraf.'5^ 


^i J .Onsager^*^  "  ^i^j(®o  +  ®2X^’  x  ^  *  (5.2) 

The  exact  trial  solution  is  of  the  form: 


c  .  (°0  +  W.  * 

S  ‘i  ^ 


(5.3) 


Symmetry  considerations  will  reduce  the  Cv%tt.'a.C  number  of 
independent  coefficients. 

Example(b2):  charged  hard  spheres  in  uniform  background 

The  Onsager-smeared  interactions  are  of  the  form:  I VM  1 


'TT  a rr-y. 


This  happens  in  general  only  for  hard  spheres  and  for  hard 
alligned  ellipsoids.  In  these  cases,  the  ASCI  solution  has  the 
same  analytic  form  as  required  of  the  PYHS  theory  at  low 
densities.  In  the  general  case,  however,  the  ASCL  solution  of 
the  PYHS  and  the  low  density  solution  may  have  different  forms, 
so  that  the  overlap-volume  basis  set  (3.7)  will  not  be  able  to 
provide  the  exact  solution  at  low  packing  fractions.  As  an 
ad-hoc  tentative  practical  solution  to  this  problem,  we  suggest 
to  employ  as  a  hard-core  basis  set  the  union,  of  the  sec  (3.7)  and 
the  overlap-volumes  of  shapes  corresponding  to  the  pa i r -e x c  1  ud e d 
regions,  i.e. 


iM] 


(5.6) 


Comment(fci):  A  method  due  to  Percus. 

Once  the  exact  analytic  form  of  the  dcf’s  is  given,  an 
alternative  to  the  variational  solution  in  provided  by  the  simple 
and  powerful  method  of  Percus^*1^.  Use  the  OZ  relations  to  expand 
the  dcf's  around  the  origin  f*«©  .  With  the  known  analytic  form 
of  the  solution  this  will  provide  the  set  of  Algebraic  equations 
for  the  coefficients.  i.e.  a  complete  analytic  solutions  without 
resort  to  the  factorization  techniques 

6.  Approximations  within  the  PY  theory  for  hard-objects. 

The  formulation  in  secs.  1,2  holds  also  for  uncharged  hard 
objects  in  the  PY  approximation.  The  main  points  to  note  is  that  T' 


H  I  i  l  >i< 


has  a  thermodynamic  seaning  1  .2  ; 


The  analogy,  pointed  out  by  Onsager,  between  a  mixture  of 
particles  of  different  orientations  and  a  mixture  of  hard  spheres 
of  diffrent  sizes,  is  born  out  manisfestly  in  our  formulation. 

In  the  asymptotic  strong  coupling  limit  we  predict  the  following 
results  for  arbitrarily  shaped  convex  hard  objects:  divergence 
at  -1  .  (3.7),  and  (3.1),  (3.2). 

Consider  now  the  following  general  approximation, 

’  C.A'i)  C.AV)  '  (6-0 


by  which  the  generating  functional  Zc  takes  the  following 
( "diagonalized” )  form: 

2  '  +  -ir  [>■*)  j<**  [*-  «  L *;  t.Ak)]  ( 

to  be  solved  via 


£  _  o  r  < 

s  C^z) 


In  the  limiting  case  when  all  objects  are  equal  size  spheres, 
eqs.  (6.2.),  (6.J)  reduce  to  the  exact  PY  equation  for  :  ard 
spheres. 


solution  to  the  Pi  equation  can  oe  judged  ay  \  s  .  g  /  cc-pa:  .0 

with  the  available  exact  solution  for  hard  sphere  mixtures,  for  a 
binary  mixture,  (6.1)  leads  to  the  following  relation  among  tne 
partial  structure  factors: 


s'  l  vw  [ 


(6.4) 


This  relation  holds  well  only  for  mixtures  of  spheres  with 
relatively  small  size  differences.  The  Pynn-Lado  21,22  ansatz 


namely 


c ^  ’  c=  (?t, 


(6.5) 


where  Cq  (  ,  r/P)  is  the  solution  of  the  PY  approximation  for  equal 
site  spheres, manifestly  violates  the  overlap-volume  analytic  form 
for  uneaual  size  spheres  (e.a.  compare  (5.1)  with  (5.5)),  and  so  cannot 
be  expected  to  be  accurate  for  mixtures  with  large  differences  in  particle 

sizes.  Lado's  numerical  results  are  encouraging  to  pursue  this 
t,ine  of  scaling  type  approximation.  It  would  have  been  of 
interest  to  compare  the  results  of  (6.5)  to  those  obtained  from  (i .  \) 
Note  that  (6.3)  is  an  exact  ASCL  result,  so  that  (<*»•  5") 
ouvncA  (."i.w^  v-e  tV\  e.  exact  results  for  the  PY  theory  near 

the  1/7  -*  \  limit: 

2^=  — — — ^  (6'c) 

7r-' 1  ‘  0  -  'Zt) 

It  is  interesting  to  observe  that  the  most  sucessful  theories  for 
the  EOS  of  mixtures  of  hard  spheres  or  non-spher ical  hard 
objects,  like  tha  ^-expansion  and  the  scaled  particle  theory, 
implicitly  assume  (6.  G  )  or  its  corresponding  "virial"  PY-EOS: 


Approximations  . 

One  of  the  aims  of  Onsager's  paper  was  to  investigate 
simplified  electrostatic  metals  of  bonding  and  aggregation.  The 
mapping  of  Onsager's  procedure  on  the  ASCL  of  integral  equation 
theories  of  fluids,  enables  us  to  take  a  further  step  in  the 
direction  pointed  out  by  Onsager. 

Since  the  Onsager  bounds  correspond  to  the  exact  ASCL  result 
of  the  HNC  and  MSA  theories,  it  is  possible  to  compare  the 
Onsager  bounds  with  the  potential  energy  of  different  possible 
structures.  To  the  extent  that  the  Onsager  energies  correspond 
to  a  unique  local  structure , they  may  predict  the  bonding  and 
aggregation  effects  within  these  integral  equation  theories  in 
strong  coupling,  without  the  need  fora  detailed  solution  of  those 
complex  equations. 

To  illustrate  this  Idea,  consider  a  binary  mixture  of  hard 
particles  of,  say,  different  sizes,  with  opposite  charges 
situated  inside  the  particles  just  off  the  surface. 

Specifically,  let  $  be  the  distance  of  each  charge  from  its 
surface  boundary,  and  let  ^ fs  •  w^ere 

are  typical  dimensions  of  the  particles.  The  Onsager  self  energy 
for  each  such  object  with  charge  Q  (or-Q)  the  surface  is 


po 
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o© 
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(7.1) 


This  is  equal,  however,  to  the  potential  energy  of  the  two 
oppositely  charged  objects  when  they  are  forming  a  molecule  , 
with  the  equal  and  opposite  charges  are  at  distance  ^  & 
from  the  other  (Fig.  Q.a).  Note  that  this  is  a  unique  local 
structure  corresponding  to  (7.1),  and  thus  we  are  able  to 
predict  that  the  HNC  and  MSA  theories  will  feature  this  bonding 
effect  when  \ 

This  simple  example  ca  be  generalized  to  local  structures 
of  higher  complexity  (e.g.  Fig.  an<*  to  include  polarization 

and  other  internal  electrostatic  effects  that  lead  to  aggregare 
formation.  The  treatment  of  polarization  in  Onsager’s  paper  can 
be  readily  incorporated  into  the  mapping  on  the  integral  equation 
theories.  It  may  be  thus  possible  to  consider  complex  structures 
as  rings,  water  etc.  Although  the  possible  aggregates  that 


>rap  lt^j 


with  Onsager's  bounds  do  not  always  represent  a  unique 


local  structure  (e.g.  charges  near  the  centers  of  the  objects) 


they  should  be  of  help  in  the  analysis  along  the  lines  set  by 
Onsager's  work,  which  -in  view  of  o\)  V  analysis-  it  is  now 
possible  to  pursue. 


8.  Further  Implications  and  Conclusions 

The  physically  intuitive  meaning  of  the  dcf's  and  its  role  in  a 
variational  solution  of  the  MSA  -HNC  type  equations  by  expansion  in  the 
"Onsager  basis  set"  ,  is  not  limited  to  uniform  systems.  It  is  a 
standard  procedure  [27],  in  the  treatment  of  a  system  of  hard  particles 
near  a  neutral  or  charged  wall  (e.g.  fluid  near  an  electrode) ,  to 
consider  a  mixture  of  different  size  particles  and  to  let  one  radius  go 
to  inf inity. This  one  particle  represents  the  hard  wall. In  such  a 
limiting  process,  done  in  the  context  of  a  specific  approximation,  such 
as  the  MSA  or  HNC  preserves  the  fundamental  form  of  the  various  dcf's 
(particle-particle, particle-wall)  retains  its  intuitive  meaning  as 
overlap  volumes  or  smeared  interactions,  and  the  procedures,  as 
discussed  in  this  paper,  can  be  carried  out  in  complete  analogy  to  the 
uniform  fluid  case. Work  along  these  lines  is  currently  under  way. 

The  Onsager  procedure  in  conjugation  with  the  Ewald  ident  ity 
serves  as  a  guideline  for  the  solution  of  the  ASCL  of  th  HNC  or  MSA 
equations  for  systems  of  concentrated  charges  without  hard  cores, 
e.g.  plasmas  ,  line  charges,  etc. (namely,  the  soft-MSA  context , [ 2 , 3 , 6 ] ) 

It  provides  the  rationale  for  the  emergence  of  the  "ion-sphere" 
boundary  conditions  in  the  treatment  of  high  density  electrons+ions 
matter  [3b]. The  insight  we  gain  from  te  general  HNC-MSA  problem  for 
HCGF  systems  may  be  used  also  in  constructing  analytic  solutions  for 
centrally  charged  spheres  for  non-uniform  systems  .Examples  of  such 
solutions  are  currently  under  investigation. 

Our  novel  method  applies  to  models  of  matter  of  a  wide  variety, 


Restricted  Primitive  Model  (RPM)  of  electrolytes. 


The  RPM  consist  of  a  mixture  of  charged  hard  spheres  of 
equal  sizes,  and  obeys  total  charge  neutrality  The 

MSA  for  this  model  has  been  solved  exactly  to  give  the  following 

Os) 

dcf  s : 


^  C  A  r  C,  C  ^  +  i.  ft  a  ^  k  +  ±  6^ 


<ncA.i) 


where  B  has  the  property 


a  _  _ 


(A. 2) 


and  CQ  (r)  is  the  PYHS  result  for  equal  size  spheres  of  diameter 

er 

In  the  limit  Qv-— >  «o  we  have 


ft 


Ct.  Cr") 


OO 


lTl*)  L 


r4q-  • 

-L  I  |  3^  <D 

1  1  '  (A. 3) 


where  ^  (.r^  is  the  electrostatic  interaction  between  two 
uniformly  surface  charged  spheres. 


Appendix  B:  Charged  hard  spheres  ift  uniform  penetrating 
background 

Defining  (b  ^'/(T  *  ^  ^  =  -E  *  S'  ,  /-/t  and  taking  the 

limit  ^  -*  for  fixed  £ <  \  ,  the  solution  of  Palmer  and 

Weeks  (PW)  takes  the  form: 


Note  that  1*) 

1.0  =» 


I  ^  » 


-f  ‘z.u+O  *  -  4  2  * 


can  be  written  as 


\Ts 


i 

q  <2  X 

(B.l) 
x  <.  \ 


(B.2) 


where 

^  (.  O  -  ^  C  y  /) 

^Sl  U>)  -  (V'?),‘  O 


(B.3) 


(B.4) 


(  B  .  5  ) 


These  functions  are,  respectively,  the  volume-volume , 
surface-surface,  and  volume-surface,  electrostatic  interactions 
between  two  spheres  of  unit  diameters,  each  composed  of  charge 


(1-rj  )  spread  uniformly  on  the  surface,  and  charge  ^  spread  uniformly 
in  the  volume. The  Onsager  object  consists  of  the  negative  of  these 
charges  plus  the  original  central  unit  point  charge. Using  (8=1)  the 
definitions  of  sec.  (4.b) 

ss-d-r^  ) 

T. 

bb=6/5  r\ 
bo=-3r^_ 
os— 2  (l-r|_  )X 

bs-2r^  (l-rj_  )  (36) 

we  find  that  ,  in  accordance  with  the  results  of  section  (4.b)  the 
internal  energy  is 

oO  J 

B(U/N)  — (1+n  -n  )8*— ss+bb+bo  (B7) 

in  agreement  with  the  result  obtained  by  direct  integration  of  (3.17). 

Appendix  C  :  Onsager  Solution  for  Hard  Spheres  with  Centrally  Imbedded 
Point  Multipoles. 


Appendix  C  :  Onsager  Solution  for  Hard  Spheres  with  Centrally  xmbedded 
Point  Multipoles. 

Consider  a  sphere  of  radius  R  with  a  centrally  imbedded  multipc 
of  order  1  .The  potential  due  to  the  central  charges  is 


(r)»  4  (t  j  q  ^J©,j)/r 

0  2«,-yvL<4~ 


$ 


(Cl) 


we 


solve  the  equation  \]  =0  inside  the  sphere  with  the  bour.da: 

condition  d  (R)«-$  (R)  (see  Sec. 4). The  solution  for  the  interior  is 

u  o 


1 1*'  C-4 


r< 


(C2) 


where 


A  — ^  R 

t*V> 


(C3 ) 


The  potential  inside  the  sphere  is 


9  (r)«  4£  2qt„YrJe'tf)  (r  "r  R  > 

lt+\  ~ 


( C4 ) 


The  normal  component  of  the  electric  field  on  the  inside  surfa 
of  the  sphere  is 


E^(R)«  4TrJ 
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4-. 


q  V  (9,*)  R 

tm  »»«  1 


:c5) 


( 2 t+l) / ( £+1)  E  (R) 


.N 


(the  cr.e  due  to  tr.e  crigma- 

density 


d  (R)  =-  Ei(R)/4ff  — o  (r) 

satisfies  the  charge  neutrality  condition 


(C6) 
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Y  (0, 

»w. 
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(C7) 


It  is  easy  to  see  that  the  result  (C5)  is  equivalent  to  a  unifo 
spread  of  the  multipole  on  the  sphere. In  the  case  of . a  sphere  with  a 
point  dipole  p  in  the  center  (pointing  in  the  direction  of  z) ,  we  hav 


o'  (R) 


*-3p  cos  0/4  T  R 
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(C8) 


The  MSA  for  the  ion-dipole  mixture  was  solved  some  years 

,23. 


ago-""’-"?.  Although  the  unequal  size  case  was  discussed  a 
complete  solution  of  this  model  is  still  not  at  hand. 

Our  present  method,  however,  allows  us  to  construct  the  DCF 
in  the  ASDL  from  purely  electrostatic  calculation  consider  first 
the  most  general  case  of  two  spheres  of  radii  b^  and  b2  and  cl\arg 
distributions  J The  center  to  center  distance  is 
.  The  interaction  energy  is 


(D.l) 


Assume  first  that  b.>b2.  Then  we  can  write  (Dl)  in  the  form 


where  cp.(.2t,£;  is  the  potential  due  Co  the  charge  smeared  on  the 
surface  of  sphere  1. 


(  D  .  3  ) 


now  for  multipolar  charge  distributions  on  the  surface  of  the 
sphere  we  have 


cJV-^o 

hTT  *»*. 


■c 


•  /*- 


(D.4) 


in  a  laboratory  fixed  reference  frame.  -  (^.'’Jis  a  vector 

pointing  to  the  direction  'o,  »  ,  &i  in  a 

reference  frame  fixed  to  molecule  i.  a^>is  the  value  of  the 
multipole  moment  of  order  m,  polarization^.  We  use  Edmonds 
notation  for. the  Wigner  spherical  harmonics 


2  S' 


w> 

£ 

•r- 


We  adopt  a  reference  frame  in  which  the  z-axis  is  the  r  vector. 


,  1-»l) " 
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/■V 

where  now  «  is  the  d  i  r  ^ f.  t.  v  o  n  i  n  3  rnoleculdr  reference  f  r  3  m 

and  the  fyUranglea  -Tl,  •*  <*.  ^  .  give  Che  orientations  o 

molecule  l  in  this  frame.  Choosing  the  origin  of  the  coordina 
at  the  center  of  1  we  vri t»  CD. 3)  in  the  form 
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with 


*  --  r 


(D.8) 


We  expand  now  the  last  factor: 


Indeed  we  will  have  two  cases 

Region  outside  1 
Region  inside  l 


we  get 
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X>  (5.10) 
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so  chat  finally 
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(5.15) 


where  the  angle  X  » is  defined 


C*n  ■»  Cr  c°"»  ^0/ 


(D.17) 


The  evaluation  of  the  integrals  (D.15)  (D.16)  is  straight 

forward.  We  quote  the  results  for  m ,  n-0,1. 

To  n-ion  interaction 

o  O  f  °  \  - 

-  G.^r'  —  [_CK -loi.)  *  lr  ( '°l-t-’oa.\-r,'1  (p.  ■.$■ 

y  «o  <3  r 
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This  form  agrees  exactly  with  the  ASCL  of  that  given  by  Hiroike  ” 

~V/hen  bys  k>x  s  la  t  then  we  get 

•  ■'  l<)  r  l 1  -  4 1  ^  C  v-  £ ']  (P  ■ ;  v 

in  agreement  with  known  result  (A3)  .  «*•  r  tbe.  Lou.  «U'^oL«. 
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and  finally  the  dipole-dipole 
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-x-*  X  <.  \ 


(5.4) 


which  is  view  of  (5.1)  is  also  the  form  of  the  exact  trial 


function . 

Example(b3):  hard  spheres  will  imbebbed  monopoles  and 

dipoles . 

For  unequal  size  spheres,  the  overlap  value  function  has  the  form 
(with  X  =  r-  >  ):l'il 


C  ^  -  A 

N — '»».  ViC. 


-  ^  V  *  *  /r  (5.5) 

The  Onsager  smeared  functions  are  given  in  Appendix  D.  when 
forming  the  trial  dcf's  we  keep  the  angular  dependence  and  attach 
coefficients  to  the  radial  part  of  each  function.  Compare  w  »-  th 
Wertheimii  solution. 

Coramen t ( b4  )  :  over  1  a p- s phe r e s  and  pa i r -e xc 1 ud ed  volumes  for 
hard  objects. 

At  low  densities,  at  the  2nd  virial  level  of  approximation,  the 
hard-core  dcf  is  given  by  the  overlap-volume  of  two  objects 
having  the  shape  of  the  pair  excluded  volume  between  the  two 
relevant  particles.  As  long  as  the  pair  excluded  volume  has  the 
same  shape  as  the  two  hard  particles,  then  the  two  kinds  of 
overlap  volumes  mentioned  above  have  the  same  analytic  form. 
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